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Abstract 

We review a formalism that includes the effects of nonperturbative 
U(l) symmetry breaking on the QCD evolution of nonsinglet struc- 
ture functions. We show that a strong scale dependence is generated 
in an intermediate energy range 0.5<^(5 2 <^5 GeV 2 for all values of x. 
We show that this explains naturally the observed violation of the 
Gottfried sum, and allows a determination of the shape of the non- 
singlet structure function, in excellent agreement with experiment. 
We argue that these effects may affect the determination of a s from 
deep-inelastic scattering. 
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Accurate measurements of nonsinglet nucleon structure functions have 
been accomplished only rather recently. Because gluons decouple from non- 
singlet quantites, such measurements provide a direct handle on the quark 
structure of the nucleon, and their scale dependence can be reliably calcu- 
lated in perturbative QCD. In particular, the first moment of the nonsin- 
glet nucleon structure function F 2 NS = i*f — F£ nas been determined rather 
accurately [JT] : 

S G = | ^ F 2 NS = 0.235 ± 0.026. (1) 

This experimental result^ is somewhat puzzling in view of the interpreta- 
tion of _F 2 NS and its first moment Sq in the QCD parton model: in the DIS 
scheme to all orders 

F 2 (x) = xY j e 2 l {q i {x)+q l {x)), (2) 

i 

where qi{x) are the distribution of quarks of flavor % in the given target, so 
that, using isospin, 

F 2 NS = I (f{x) + *P(x)) - (<F(x) + <?(x))) , (3) 

where u p and d p are up and down quark distributions in the proton. Now, 
Sq is to leading order scale independent; hence at all scales Sq would be 
expected to be close to its quark model value, obtained assuming the quark 
distributions in Eq. (^) to be given by valence quarks: this, however, leads to 
S G = |. 

A possible interpretation of this discrepancy is that isospin is broken; 
however, isospin violation in this channel is Q at least one order of magni- 
tude smaller than required to explain Eq. ([[]). Thus the data (if correct) imply 
that it is the identification of the quark distributions in Eq.(|3]) with valence 
distributions to fail: the sea must have a large nonsinglet component, anticor- 
related to the valence. It is easy to reproduce this sort of scenario in effective 
models of the nucleonQ. However, due to the (one-loop) scale independence 
of Sg, such a nonsinglet component cannot be generated perturbatively from 
the starting valence distributions. 



*Very recent data[^) show evidence for shadowing in scattering on deuterium. If this 
effect is confirmed the value Eq.(|l]) should be reduced by a further 10-15 % 
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The evolution equation for nonsinglet quark distributions takes the simple 
general form|| 

j t [q±q) = (Q gq ±Q qi} )®[q±q), (4) 
to any perturbative order, where 

q = u-d- Q qq = vl-Vl D , (5) 

and V D (V ND ) is any quark-quark splitting function V m such that i = j 
(i 7^ j). At one loop, the only process which leads to nonsinglet evolution 
is gluon radiation (Fig. la), which contributes to P qq but has vanishing first 
moment. This conservation holds true to all orders in the charge-conjugation 
odd case, but fails beyond leading order in the C-even one. Two-loop evolu- 
tion is then driven by the diagram of Fig. lb: whereas the flavor-invariance 
of QCD leads one to expect that this should contribute equally to V D and 
, p ND ; due to Fermi statistics the final state must be antisymmetrized with 
respect to the two identical quarks when i = j, while this is not the case if 
i 7^ j. This is enough || to lead to nonvanishing values for all the nonsinglet 
anomalous dimensions, which, however, remain very small. 

In QCD, axial U(N/) flavor symmetry is broken down to SU(N/) by non- 
perturbative effects due to the axial anomaly. This leads to a large difference 
between flavor-preserving and flavor-changing transitions 0], and thus to siz- 
able values of the anomalous dimensions |7]]. Indeed, assume that the emitted 
and unobserved q-q pair may form a bound state (Fig.lc). Emission of 
neutral (charged) mesons will then contribute to V qq (V qq D ). In the pseu- 
doscalar sector, where U(l) symmetry breaking manifests itself, 7r ± emission 
will thus contribute to V qq D , and tt° and r\ emission to V qq (considering for 
simplicity the SU(2) case). Because of the much larger mass of the i] the 
flavor-preserving process is suppressed, thus leading to a large negative value 
for Q qq . This leads to sizable evolution, with the sign required to explain the 
result Eq.(p. 

This argument can be made quantitative |7|] by setting up generalized 
evolution equations which include an effective coupling to bound states of 
the form of Fig.lc. Because symmetry breaking only appears in the Gold- 
stone sector, only pseudoscalar mesons need to be included. Introducing a 
nonsinglet pion distribution tt = 7r + — 7r~, the nonsinglet evolution equations 
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are then 



-^q = Qqq®q+ Qgg ® <? + V q7T <g> 7T, (6) 

d _ 

^g = Qgg <S> g + Qgg <E> g - V q7T <g> tt, (7) 
= P vq <B,(q-q)+V mr ®n. (8) 

It is still necessary to specify the effective quark-meson coupling, denoted 
by a blob in Fig.lc. The total cross section for the process of Fig.lc turns 
out to depend only on a pseudoscalar and an axial vertex function; the 
former dominates for intermediate values of Q 2 (up to a few GeV 2 ), while 
the latter controls the large Q 2 tail. As a a consequence, the bulk of the 
splitting function depends on one single parameter A, the radius of the form 
factor, which is related by the chiral Ward identity to the constituent quark 
mass M q . The large Q 2 tail depends also on the axial coupling g n , which is 
expected to be 

It may now be checked explicitly that V qq can be computed in the LLA, 
in that the cross section is dominated by a collinear singularity where the 
coupling is effectively pointlike. The anomalous dimensions are sizable for 
values of Q 2 up to 5-10 GeV 2 , while at both very small Q 2 ~ 0.05 GeV 2 , 
and at large Q 2 > 10 GeV 2 , all the anomalous dimensions flatten, ensuring 
smooth connection to a valence quark picture in the infrared and to the 
usual perturbative behaviour in the ultraviolet. As expected, , p Nr> is larger 
than 'p Nr> by roughly a factor two , thus leading to significant nonsinglet 
evolution. 

The scale dependence of Sq can now be determined: the result obtained 
assuming Sq to take the "quark model" value Sq = 1/3 at a reference scale 
Qo ~ 200 MeV is shown in Fig. 2 as a function of the two parameters M q 
and <jv, and compared to the experimental result. Very good agreement is 
obtained with reasonable values of the parameters. This means that the full 
nonsinglet sea is generated dynamically: therefore, not only its first moment 
but its full shape can now be predicted 0. 

To this purpose, we need a model for the starting valence and its evolution 
at low scale. Once this is specified 0, the nonsinglet structure function is fully 
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determined at all x and Q 2 , because the bulk of the low-Q 2 nonperturbative 
effects is expected to be encapsulated by the anomalous evolution equations 
(0~§P- The scale dependence is driven by two competing mechanisms: gluon 
radiation (Fig. la) which softens the valence quark distribution, and bound 
state emission (Fig.lc) which leads to to a net increase in the number of 
sea q-q pairs. Thus, the peak of the nonsinglet distribution as a function 
of x remains approximately fixed at its starting value x ~ |, rather than 
shifting towards smaller x as it happens in the singlet case. However, the 
overall decrease of the structure function is stronger than found with ordinary 
perturbative evolution. 

The structure function determined thus is compared to the dataQ and to 
the result of ordinary perturbative evolution in Fig.3. It is apparent that the 
effects of anomalous evolution are large for all values of x. The (unphysical) 
effect of switching gluon emission off is also shown. In the large-x region, 
where the NMC data have large errors, an alternative comparison can be 



made with more precise data at a larger value of Q Oj. The stability of 



these results upon variation of the parameters has been checked explicitly 



Drell-Yan scattering can provide a direct determination of 

= d{x) - ujx) 
l J_ d(x)+u(x)- 1 ] 

Predicting the value of R(x) in the anomalous evolution scenario would re- 
quire a computation of singlet evolution as well; we may however estimate 
it by assuming the bulk of the effect of anomalous evolution on the singlet 
to be just the production of an asymmetric sea. The result (Fig. 4) has the 
generic feature of saturating to a constant when x^0.2 because at large x 
the symmetric sea vanishes more rapidly than the anomalously generated 
asymmetric one (as discussed above), thus only the latter contribution to R 
survives, which is just equal to a constant group-theoretic factor. Our result 



is in excellent agreement with the data |fl2| , |T3f| , which however cannot ex- 
clude some of the alternative models; a measurement of R Eq. (||) for several 
values of x would be required to discriminate between them. 



^Notice that the data have been recently revised; the new data agree better with our 
prediction [|| than those pj to which we originally H compared it. 
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The scale dependence of the full nucleon structure function is dom- 
inated by nonsinglet evolution if x^0.3 because the gluon distribution falls 
very rapidly at large x. Thus, anomalous evolution will affect the scale de- 
pendence of F$ at large x and intermediate values of Q 2 . We may estimate 
this effect by again assuming that the bulk of the anomalous contribution to 
both singlet and nonsinglet is just due to the diagram of Fig.lc. Then, in 
this range of x and Q 2 Ff will decrease less rapidly as a function of log Q 2 
then it would if anomalous evolution were neglected, essentially because of 
the anomalous production of q-q pairs discusse above, which counters the 
softening of the structure function (Fig. 5a). Of course, at large Q 2 the usual 
results are regained. 

A decrease of \dlnF2 /d\nQ 2 \ of the same size is found if a s (M z ) is re- 
duced by roughly 10 % (Fig5b, from Ref.(T^|). Because the value of ot s is 



extracted |T5|] from F 2 by fitting its scale dependence while neglecting anoma- 
lous evolution, if the effect is indeed present the value of a s is systematically 
underestimated. Of course the relevant fits include also the small-Q 2 re- 
gion, where higher- twist effects are dominant, and the large-Q 2 region, were 
anomalous evolution is negligible; a new global fit is thus required to pin 
down the size of the effect. The sign of the effect and order of magnitude, 
however, appear to be just what required in order to bring this determination 
of a s in line with the LEP value. 

Acknowledgements: Part of the work described here was done in col- 
laboration with V. Barone and M. Genovese. We thank G. Altarelli for 
pointing out the relevance of our results to the determination of a s , and 
M. Arneodo and C. Quigg for discussions. 
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Figure Captions 



Fig. [1] Diagrams which contribute to V qq : a) perturbative one loop; b) 
perturbative two loops; c) one loop bound-state emission 

Fig. [2] Anomalous evolution of Sq (Eq.0) for various values of the con- 
stituent quark mass M q (a) and of the axial coupling g n (b). 

Fig. [3] Comparison of the nonsinglet structure function generated dy- 
namically by gluon and bound state emission with the data|lj]. The re- 
sults obtained switching off gluon emission or bound state emission are 
also shown. 

Fig. [4] Comparison of the experimental results for the SU(2) antiquark 
asymmetry R (Eq.§) with the prediction of various models. 

Fig. [5] Change in the scale dependence of the structure function F2 due 
to the inclusion of nonsinglet anomalous evolution effects (a) or variation 
of the value of a s (b, from Ref.||T5|). 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409249v3 
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Fig. 1 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409249v3 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409249v3 
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This figure "figl-4.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409249v3 
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This figure "figl-5.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409249v3 
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